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ABSTRACT 

The synchronizing of a wind turbine generator against an infinite 
bus under random conditions is studied for the first time. With a digi- 
tal computer, complete solutions for rotor speed, generator power angle, 
electromagnetic torque, wind turbine torque, wind turbine blade pitch 
angle, and armature current are obtained and presented by graphs. Exper- 
iments have been recently performed on the ERDA-NASA 100 kW wind turbine. 
Experimental results matched computer study results very closely and con- 
firmed that the synchronization can be accomplished by means of the exist- 
lng — sj^eeh^ontf ol“ "sys tem~'and~ ah automatic synchronizer . 


INTRODUCTION 

The ERDA-NASA 100 kW wind turbine is the first phase of the ERDA wind 
energy program performed at NASA-Lewis Research Center. This turbine, 
designated the Mod-0, became operational in September, 1975, at the NASA 
Plum Brook Station near Sandusky, Ohio. This project and its present 
status are described in a recent publication [1] . 

The Mod-0 wind turbine (Fig. 1) is a horizontal-axis, propeller-type 
machine. A two-bladed, 125-foot diameter rotor, mounted on a 100-foot 
truss tower, drives an alternator through a step-up gear box. The rotor 
is downwind of the tower and rotates at a constant speed of 40 rpm. The 
alternator is a 125 kVA, 0.8 powev factor, three-phase, Y-connected syn- 
chronous machine. Figure 2 shows details of the wind turbine system and 
assembly. 

t- 

One of the objectives of the Mod-0 project is to gain early experi- 
ence with wind turbines operating in utility networks [1]. Therefore, 
i synchronization becomes an important subject. However, the synchronizing 
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of a wind turbine with a large utility network has not been fully studied. 
The purpose of the work reported here is to investigate the system 
dynamics following the synchronization and to determine the feasibility 
and optimum conditions for achieving the required operation. 

A mathematical model is developed first, according to the specifi- 
cations and the configuration of the wind turbine system. Then, by using 
a digital computer and the mathematical model, studies are performed to 
evaluate the variations of rotor speed, turbine torque, and transients in 
the alternator following the synchronization. The initial frequency and 
phase position of the alternator's voltage are assumed to be different 
from those of the bus voltage. Wind gusts are neglected. Later, it is 
intended to study the effect of strong wind gusts on the performance of 
synchronization. However, the simulation model presented here can be 
readily expanded to include the effects of strong wind gusts. Recently, 
experiments were performed on the Mod-0 wind turbine under fairly steady 
wind speed condition. Both experimental and computer study results con- 
firmed that synchronization of a wind turbine can be accomplished if the 
properly designed speed control system and synchronizer are used. The 
developed digital computer program for simulating wind turbine synchroni- 
zation should prove useful for future related studies as well. 


SYSTEM EQUATIONS 

Figure 3 is a schematic of the system under study. Details of a 
representative speed control system are shown in Fig. 4. The alternator, 
the mechanical dynamics, and the speed control system are represented by 
equations given in the following section. 


Synchronous Machine Equations 

The alternator is the most complex component of the generating system 
and largely dominates the performance of synchronizing. A fairly rigorous 
representation is required. 

Park's equations [2] can be used in the analysis, but they are not 
given in a form suitable for solution by means of digital computers. 

Also, some constants used in Park's equations cannot be measured by 
simple tests. Park's equations were modified by Olive [3], and rearranged 
in this paper. The following equations can be easily verified, when nec- 
essary. 
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The armature current is given by Park [2] : 


i » i , cos tot - i sin tot (11) 
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The electromagnetic torque opposite to the direction of rotation is : 
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Equations (1) to (12) are in per unit based upon the machine rating. 
These equations will include the effect of external series impedance, if 
the machine impedances are replaced by the sum of the external impedance 
and machine impedance. 
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Equations of Motion 

The electromagnetic torque (referred, to the rotor shaft) in lb-ft is 


T t = (kV base) (738)T 

L L to e 

The equations of motion, or swing equations, are 
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where S m is in mechanical radians and T in lb-ft. 

The torque angle of the alternator in electrical degrees is given by 


6 = 180 NP 6 
tt 2 m 


(16) 


is 


The angular speed of the alternator in electrical radians per second 
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Equations of the Speed Control System 

The effects of the speed governor control during synchronizing per- 
iods can be taken into account by using the simplified representation of 
the speed control system shown in Fig. 4. This representation includes 
transfer functions describing the pitch change mechanism and a transfer 
function describing the blade torsional dynamics. The input and output 
variables of these transfer functions are indicated in Fig, 4. The input 
and output variables of the flow limiter are: 

X = (X n - K,n - X.) (18) 

T i 1 e L 
P 


Y - limit {^- 22 _ 5 , 22 5 , xj 


(19) 


Other six equations can be easily formulated and are not given here. 
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The corresponding torque developed by the rotor in Ib-ft is 

t = £(v w ,n rI 0) ( 20 ) 

Equation (20) is a nonlinear function of wind speed in mph, rotor 
speed in rpra, and pitch angle in mechanical degrees [4]. The turbine 
torque-pitch angle curves are given in Fig. 5. The blade pitch angle 
physically can vary only from -90 degrees to zero degree sc that the 
slope of the torque-angle curves is always positive. As one per unit 
torque referred to the rotor shaft of the Mod-0 wind turbine is equal 
to 22 023 foot-pounds, the maximum turbine torque, at a 15 mph wind speed, 
is less than 15 000 foot-pounds or less than 0,7 per unit. 


DIGITAL SIMULATION 

In this study, there are totally 26 nonlinear system equations. 

They are Interrelated, and therefore must be solved simultaneously. An 
analytical solution is impossible, but a digital simulation will give the 
required information. The fourth -order Adams -Bashforth (predictor) and 
Adams -Moulton (corrector) method is used as a mathematical tool to obtain 
the required solutions. This method is popularly used, because it is 
less laborious than other methods and yet gives comparable accurate 
results. The predictor-corrector method in this study needs only 1/3 
of the computer time required by the fourth-order Runge-Kutta method for 
the same simulation. 

In terms of d, q components, the initial voltages across the 
synchronizing breaker contacts are 

V, = sin 5 — sin (to - to)t (21) 

do o 

V = -cos 6 + cos (a) - co)t (22) 

qo o 

The effect of closing the breaker contacts is simulated by applying 
— V{j 0 and -V qo in (1) and (2), respectively, and setting the field 
voltage Vf equal to zero in (3). The initial value of field current is 
included in the numerical computations of electromagnetic torque given by 
( 12 ) . 

In the analysis, the mechanical and aerodynamic losses are neglected 
in order to study the problem under worst possible conditions. The effect 
of neglecting these losses may be to increase the magnitude and duration 
of transients during the synchronization process. 
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System Constants 
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The initial pitch angle 6, corresponding to the wind speed, initial 
turbine torque, and initial rotor speed is calculated by using the torque 
routine [4]. For each speed error the Initial generator power angle or 
phase position of the alternator relative to the infinite bus voltage is 
calculated according to the specifications of the synchronizer and cir- 
cuit breakers used in the Mod-0 generator system. Other initial condi- 
tions can be readily determined. 


RESULTS OF COMPUTER STUDIES 

In this study, the two extreme upper speed limits of closing the 
existing synchronizer contacts are used. They are 7.5 and 2.5 rpm above 
1800 rpm (synchronous speed of the synchronous machine). 

Procedures of bringing the machine to these speeds before closing the 
synchronizing switch are outlined in a previous publication [5], A 
15 mph wind speed and the 0.02 per unit reactance of the transformer and 
the tie line (Fig. 3) are also used in the simulation. The following 
results are observed: 

Figure 6(a). - The initial rotor speeds referred to the low speed 
shaft are 40.1667 and 40.0555 rpm, respectively. For thfc first case 
(initial speed error =7.5 rpm on 1800 rpm basis), the rotor speed 
reaches synchronous speed (40 rpm) with an error less than 0.1 percent 
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within the first three seconds following the synchronization switching. 

For the second case (initial speed error = 2.5 rpm) , the correspond- 
ing settling time is less than 2 seconds. Practically, the rotor speed 
for either case will reach synchronous speed completely within first 2 or 
3 seconds after closing the synchronizing switch as confirmed by experi- 
ments. The ripples of the rotor speed curves would completely vanish 
within few seconds if the system damping is included in the simulation. 

Figure 6(b) . - The plot of generator power or torque angle versus 
time or the swing curve is the best indication of the system stability 
under transient conditions. The swing curve for either case shows that 
a system equilibrium, or a smooth synchronization will ultimately be 
accomplished. For the first case, the peak valu? of the power angle 
amounts to only about 27 electrical degrees (t = 0.4 sec). The power 
angle, at the end of the simulation, is less than 2.2 electrical degrees. 
For the second case, the peak and final values of the power angle are 
about 15 electrical degrees (at t = 0.3 sec) and about 1.4 electrical 
degrees (at t = 10 sec) , respectively. The small variation of the power 
angle in each case assures that the operation of synchronization under 
assumed conditions is satisfactory. 

Figure 6(c) . - For the first case, the maximum and minimum values of 
electromagnetic torque are 2898 ft-lbs (at t = 1.28 sec) and -3076 ft-lbs 
(at t = 3.02), respectively. They are about 0.13 per unit and -0.14 per 
unit base upon the machine rating. The positive torque is a generator 
torque and the negative torque is a motor torque. For the second case, 
the corresponding values are equal to about O.i per unit (at 
t = 1.40 sec) and -0.1 per unit (at t = 2.92 sec) , respectively. In 
each case, the magnitude of the electromagnetic torque is reasonably 
small. This fact further demonstrates that the synchronizing can be 
completed smoothly. 

Figure 6(d) . - The variation of wind turbine torque in each case is 
consistent with the changes of the electromagnetic torque and the rotor 
speed. At any instant, the three torques, namely, electromagnetic torque, 
wind turbine torque, and inertia torque must satisfy the equation of 
motion (15). The negative wind turbine torques are caused by the high 
rate of decreasing the rotor speeds at the corresponding moments. 

Figure 6 (e) . - The curve given in this figure shows the variation of 
wind turbine blade pitch angle. The wind turbine torque and the rotor 
speed are controlled by varying the blade pitch angle. 

Finally, the armature currents are calculated for both cases. How- 
ever, the curves are not given here as they cannot be reproduced clearly. 
For the first case, the maximum and minimum values of the armature cur- 
rents are 0.785 per unit (at t = 2.26 sec) and -0.816 per unit (at 
t = 2.30), respectively. At the end of the simulation, the current if, 
equal to -0.223 per unit. For the second 


8 


case, the corresponding values are 0.607 per unit (at t “ 3.3 sec), 
-0.627 per unit (at t = 3.24 sec), and -0.031 per unit (at t = 10 sec). 
The small currents indicate that the transients caused by the synchroni- 
zation switching are limited. 


CONCLUSIONS 

1. For both selected initial rotor speeds, experimental and computer 
study results consistently confirmed that the required synchronization 
can be accomplished with limited transients by means of an effective speed 
control system and an automatic synchronizer. The initial rotor speed 
closest to synchronous speed aroears to be more favorable for a smooth 
synchronization, as the corres' ndlng transient is relatively smaller, 

2. The computer program developed for the analysis can be used as a 
basis for determining specifications and design criteria for future wind 
turbine synchronization systems. The program results are conservative 
in that the synchronization has been studied for the worst possible con- 
ditions . 

3. Wind energy is unsteady, particularly at the low wind speeds at 
which the wind turbine will generally be brought on line. The computer 
program can be readily modified to study synchronization phenomena in 
gusting wind conditions. 


NOMENCLATURE 


Synchronous Machine Constants 
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normal frequency = 60 Hz 

direct axis and quadrature axis currents 

field current 

self-inductance of field circuit 

self-inductances of direct- and quadrature-axis damper 
circuits 

mutual inductance between direct-axis armature circuit and 
field circuit 

mutual Inductance between direct-axis armature circuit and 
direct-axis (or quadrature-axis) damper circuit 
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armature resistance per phase 
machine electromagnetic torque 

open circuit field sub transient and transient time con- 
stants 

open circuit subtransient time constant i>f quadrature- 
axiB damper circuit 

terminal voltage of Infinite bus 

applied voltage to field circuit 

applied voltages to direct-axis and quadrature-axis arma- 
ture circuits , respectively 

limiting reactance including the 0.02 p.u. transformer 
reactance 

synchronous reactance in direct and quadrature axis cir- 
cuits, respectively 

direct-axis transient reactance 

direct-axis subtransient reactance 

torque (or power) angle in electrical radians 

flux linkages of direct-axis and quadrature-axis armature 
circuits , respectively 

flux linkages of field, direct-axis, and quadrature-axis 
damper circuits 

instantaneous angular speed in electrical radlans/sec 
base angular speed in electrical radlans/sec = 2rf 


Mechanical Constants 
B 


damping coefficient of the complete rotating system in 
lb-ft-sec 

inertia of the complete rotating system (referred to low 
speed shaft) in lb-ft-sec^ 
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gear ratio between high and low speed shaft 
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Speed 

K 1 

k 2 
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T 

E 
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torque developed by wind turbine rotor in lb-ft 
wind speed in mph 

angular displacement of wind turbine rotor in mechanical radians 

angular speed of wind turbine rotor in radians per second and in 
rpm, respectively 

initial generator rotor speed error (rpm) on 1800 rpm basis 

reference speed (referred to wind turbine rotor shaft) in radians 
per second 

Control System 

controller proportional gain in seconds 
controller integral gain per second 
pitch angle of blade in mechanical radians 
reference pitch angle of blade in mechanical radians 
damping ratio 

control hydraulic actuator time constant in seconds 
rotor rotational speed error in radians per second 
undamped natural angular frequency in radians per second 
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Figure 1. - ERDA* NASA 100-kW wind turbine. 
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Figure 2. - The 100-kW wind turbine system and as- 
sembly. 
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Figure 3. - A schematic e: the system under study. 



Figure 4. - The speed control system of ERDA-NASA 100-kW wind turbine. 






Figure 5. - The turbine torque and pitch angle curves. 
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Figure 6. - Concluded. 
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